Enhanced full-inversion-based ultrasound elastography for evaluating tumor response to neoadjuvant chemotherapy in patients with locally advanced breast cancer.

PURPOSE: An enhanced ultrasound elastography technique is proposed for early assessment of locally advanced breast cancer (LABC) response to neoadjuvant chemotherapy (NAC). METHODS: The proposed elastography technique inputs ultrasound radiofrequency data obtained through tissue quasi-static stimulation and adapts a strain refinement algorithm formulated based on fundamental principles of continuum mechanics, coupled with an iterative inverse finite element method to reconstruct the breast Young's modulus (E) images. The technique was explored for therapy response assessment using data acquired from 25 LABC patients before and at weeks 1, 2, and 4 after the NAC initiation (100 scans). The E ratio of tumor to the surrounding tissue was calculated at different scans and compared to the baseline for each patient. Patients' response to NAC was determined many months later using standard clinical and histopathological criteria. RESULTS: Reconstructed E ratio changes obtained as early as one week after the NAC onset demonstrate very good separation between the two cohorts of responders and non-responders to NAC. Statistically significant differences were observed in the E ratio changes between the two patient cohorts at weeks 1 to 4 after treatment (p-value < 0.001; statistical power greater than 97%). A significant difference in axial strain ratio changes was observed only at week 4 (p-value = 0.01; statistical power = 76%). No significant difference was observed in tumor size changes at weeks 1, 2 or 4. CONCLUSION: The proposed elastography technique demonstrates a high potential for chemotherapy response monitoring in LABC patients and superior performance compared to strain imaging.

A novel tissue mechanics-based method for improved motion tracking in quasi-static ultrasound elastography.

PURPOSE: Most cancers are associated with biological and structural changes that lead to tissue stiffening. Therefore, imaging tissue stiffness using quasi-static ultrasound elastography (USE) can potentially be effective in cancer diagnosis. USE techniques developed for stiffness image reconstruction use noisy displacement data to obtain the stiffness images. In this study, we propose a technique to substantially improve the accuracy of the displacement data computed through ultrasound tissue motion tracking techniques, especially in the lateral direction. METHODS: The proposed technique uses mathematical constraints derived from fundamental tissue mechanics principles to regularize displacement and strain fields obtained using Global Ultrasound Elastography (GLUE) and Second-Order Ultrasound Elastography (SOUL) methods. The principles include a novel technique to enforce (1) tissue incompressibility using 3D Boussinesq model and (2) deformation compatibility using the compatibility differential equation. The technique was validated thoroughly using metrics pertaining to Signal-to-Noise-Ratio (SNR), Contrast-to-Noise-Ratio (CNR) and Normalized Cross Correlation (NCC) for four tissue-mimicking phantom models and two clinical breast ultrasound elastography cases. RESULTS: The results show substantial improvement in the displacement and strain images generated using the proposed technique. The tissue-mimicking phantom study results indicate that the proposed method is superior in improving image quality compared to the GLUE and SOUL techniques as it shows an average axial strain SNR and CNR improvement of 44% and 63%, and lateral strain SNR and CNR improvement of 130% and 435%, respectively. The results of the phantom study also indicate higher accuracy of displacement images obtained using the proposed technique, including improvement ranges of 7-84% and 26-140% for axial and lateral displacement images, respectively. For the clinical cases, the results indicate average improvement of 48% and 64% in SNR and CNR, respectively, in the axial strain images, and average improvement of 40% and 41% in SNR and CNR, respectively, in the lateral strain images. CONCLUSION: The proposed method is very effective in producing improved estimate of tissue displacement and strain images, especially with the lateral displacement and strain where the improvement is highly remarkable. While the method shows promise for clinical applications, further investigation is necessary for rigorous assessment of the method's performance in the clinic.

An Enhanced Method for Full-Inversion-Based Ultrasound Elastography of the Liver.

Similar to many other types of cancer, liver cancer is associated with biological changes that lead to tissue stiffening. An effective imaging technique that can be used for liver cancer detection through visualizing tissue stiffness is ultrasound elastography. In this paper, we show the effectiveness of an enhanced method of quasi-static ultrasound elastography for liver cancer assessment. The method utilizes initial estimates of axial and lateral displacement fields obtained using conventional time delay estimation (TDE) methods in conjunction with a recently proposed strain refinement algorithm to generate enhanced versions of the axial and lateral strain images. Another primary objective of this work is to investigate the sensitivity of the proposed method to the quality of these initial displacement estimates. The strain refinement algorithm is founded on the tissue mechanics principles of incompressibility and strain compatibility. Tissue strain images can serve as input for full-inversion-based elasticity image reconstruction algorithm. In this work, we use strain images generated by the proposed method with an iterative elasticity reconstruction algorithm. Ultrasound RF data collected from a tissue-mimicking phantom and in-vivo data of a liver cancer patient were used to evaluate the proposed method. Results show that while there is some sensitivity to the displacement field initial estimates, overall, the proposed method is robust to the quality of the initial estimates. Clinical Relevance- Improved elasticity images of the liver can aid in achieving more reliable diagnosis of liver cancer.

Effect of Cardiac Phase on Cardiac Output Index Derived from Dynamic CT Myocardial Perfusion Imaging.

Purpose: The aortic time-enhancement curve obtained from dynamic CT myocardial perfusion imaging can be used to derive the cardiac output (CO) index based on the indicator dilution principle. The objective of this study was to investigate the effect of cardiac phase at which CT myocardial perfusion imaging is triggered on the CO index measurement with this approach. Methods: Electrocardiogram (ECG) gated myocardial perfusion imaging was performed on farm pigs with consecutive cardiac axial scans using a large-coverage CT scanner (Revolution, GE Healthcare) after intravenous contrast administration. Multiple sets of dynamic contrast-enhanced (DCE) cardiac images were reconstructed retrospectively from 30% to 80% R-R intervals with a 5% phase increment. The time-enhancement curve sampled from above the aortic orifice in each DCE image set was fitted with a modified gamma variate function (MGVF). The fitted curve was then normalized to the baseline data point unaffected by the streak artifact emanating from the contrast solution in the right heart chamber. The Stewart−Hamilton equation was used to calculate the CO index based on the integral of the fitted normalized aortic curve, and the results were compared among different cardiac phases. Results: The aortic time-enhancement curves sampled at different cardiac phases were different from each other, especially in the baseline portion of the curve where the effect of streak artifact was prominent. After properly normalizing and denoising with a MGVF, the integrals of the aortic curve were minimally different among cardiac phases (0.228 ± 0.001 Hounsfield Unit × second). The corresponding mean CO index was 4.031 ± 0.028 L/min. There were no statistical differences in either the integral of the aortic curve or CO index among different cardiac phases (p > 0.05 for all phases).

Measurement of the hyperelastic properties of 72 normal homogeneous and heterogeneous ex vivo breast tissue samples.

The mechanical properties of normal soft tissues, including breast tissue, have been of interest to the biomedical research community as there are many clinical and industrial applications that can benefit from quantitative information characterizing such properties. For instance, computer assisted surgery planning, elastography for breast cancer diagnosis, and bra design can all involve biomechanical modeling of the breast to predict its deformation or stress distribution. It is known that most biological soft tissues, including breast tissue, exhibit nonlinear mechanical response over large strains. As such, it is necessary to model such tissues as hyperelastic. In this work, we used indentation testing to estimate the hyperelastic parameters of 4 models (3rd order Ogden, 5-term polynomial, Veronda-Westman and Yeoh) estimated from 72 healthy ex vivo breast tissue samples covering adipose, fibroglandular, and mixed tissue. All estimated parameter sets were confirmed to represent stable material using Drucker's stability criterion. We observed that all three tissue types were statistically similar solidifying the use of homogenous breast modelling over large strain simulation.

Combining First- and Second-Order Continuity Constraints in Ultrasound Elastography.

Ultrasound elastography is a prominent noninvasive medical imaging technique that estimates tissue elastic properties to detect abnormalities in an organ. A common approximation to tissue elastic modulus is tissue strain induced after mechanical stimulation. To compute tissue strain, ultrasound radio frequency (RF) data can be processed using energy-based algorithms. These algorithms suffer from ill-posedness to tackle. A continuity constraint along with the data amplitude similarity is imposed to obtain a unique solution to the time-delay estimation (TDE) problem. Existing energy-based methods exploit the first-order spatial derivative of the displacement field to construct a regularizer. This first-order regularization scheme alone is not fully consistent with the mechanics of tissue deformation while perturbed with an external force. As a consequence, state-of-the-art techniques suffer from two crucial drawbacks. First, the strain map is not sufficiently smooth in uniform tissue regions. Second, the edges of the hard or soft inclusions are not well-defined in the image. Herein, we address these issues by formulating a novel regularizer taking both first- and second-order derivatives of the displacement field into account. The second-order constraint, which is the principal novelty of this work, contributes both to background continuity and edge sharpness by suppressing spurious noisy edges and enhancing strong boundaries. We name the proposed technique: Second-Order Ultrasound eLastography (SOUL). Comparative assessment of qualitative and quantitative results shows that SOUL substantially outperforms three recently developed TDE algorithms called Hybrid, GLUE, and MPWC-Net++. SOUL yields 27.72%, 62.56%, and 81.37% improvements of the signal-to-noise ratio (SNR) and 72.35%, 54.03%, and 65.17% improvements of the contrast-to-noise ratio (CNR) over GLUE with data pertaining to simulation, phantom, and in vivo tissue, respectively. The SOUL code can be downloaded from code.sonography.ai.

In-vivo lung biomechanical modeling for effective tumor motion tracking in external beam radiation therapy.

Lung cancer is the most common cause of cancer-related death in both men and women. Radiation therapy is widely used for lung cancer treatment; however, respiratory motion presents challenges that can compromise the accuracy and/or effectiveness of radiation treatment. Respiratory motion compensation using biomechanical modeling is a common approach used to address this challenge. This study focuses on the development and validation of a lung biomechanical model that can accurately estimate the motion and deformation of lung tumor. Towards this goal, treatment planning 4D-CT images of lung cancer patients were processed to develop patient-specific finite element (FE) models of the lung to predict the patients' tumor motion/deformation. The tumor motion/deformation was modeled for a full respiration cycle, as captured by the 4D-CT scans. Parameters driving the lung and tumor deformation model were found through an inverse problem formulation. The CT datasets pertaining to the inhalation phases of respiration were used for validating the model's accuracy. The volumetric Dice similarity coefficient between the actual and simulated gross tumor volumes (GTVs) of the patients calculated across respiration phases was found to range between 0.80 ± 0.03 and 0.92 ± 0.01. The average error in estimating tumor's center of mass calculated across respiration phases ranged between 0.50 ± 0.10 (mm) and 1.04 ± 0.57 (mm), indicating a reasonably good accuracy of the proposed model. The proposed model demonstrates favorable accuracy for estimating the lung tumor motion/deformation, and therefore can potentially be used in radiation therapy applications for respiratory motion compensation.

Characterizing regional myofiber damage post acute myocardial infarction using global optimization.

Medical imaging derived cardiac biomechanical models offer a wealth of new information to be used in diagnosis and prognosis of cardiovascular disease. A noteworthy feature of such models is the ability to predict myofiber contraction stresses during acute or chronic ischemic events. Current techniques for heterogeneous contraction models require tissue motion tracking capabilities which are neither available on all imaging modalities, nor currently used in the clinic. Proposed in this article is a proof of concept of a tissue tracking independent technique focused on shape optimization to predict the contraction stresses of in-silico left ventricle models simulating various acute myocardial infarction events. The technique involves three variables defined in the left ventricle muscle. Two of the variables represent the contraction stresses in the healthy and infarct regions while the third is a novel periinfarct variable defining a non-contracting myofiber state allowing finer classification of local myofiber damage. Results indicate that the contraction stress reconstruction errors are overall smaller than 12% when considering standard errors associated with population modelling for the new variable of interest.

MR and ultrasound cardiac image dynamic visualization and synchronization over Internet for distributed heart function diagnosis.

Dual-modality 4D cardiac data visualization can convey a significant amount of complementary image information from various sources into a single and meaningful display. Even though there are existing publications on combining multiple medical images into a unique representation, there has been no work on rendering a series of cardiac image sequences, acquired from multiple sources, using web browsers and synchronizing the result over the Internet in real time. The ability to display multi-modality beating heart images using Web-based technology is hampered by the lack of efficient algorithms for fusing and visualizing constantly updated multi-source images and streaming the rendering results using internet protocols. To address this practical issue, in this paper we introduce a new Internet-based algorithm and a software platform running on a Node.js server, where a series of registered cardiac images from both magnetic resonance (MR) and ultrasound are employed to display dynamic fused cardiac structures in web browsers. Taking advantage of the bidirectional WebSocket protocol and WebGL-based graphics acceleration, internal cardiac structures are dynamically displayed, and the results of rendering and data exploration are synchronized among all the connected client computers. The presented research and software have the potential to provide clinicians with comprehensive information and intuitive feedback relating to cardiac behavior and anatomy and could impact areas such as distributed diagnosis of cardiac function and collaborative treatment planning for various heart diseases.

Superviscous properties of the in vivo brain at large scales.

There is growing awareness that brain mechanical properties are important for neural development and health. However, published values of brain stiffness differ by orders of magnitude between static measurements and in vivo magnetic resonance elastography (MRE), which covers a dynamic range over several frequency decades. We here show that there is no fundamental disparity between static mechanical tests and in vivo MRE when considering large-scale properties, which encompass the entire brain including fluid filled compartments. Using gradient echo real-time MRE, we investigated the viscoelastic dispersion of the human brain in, so far, unexplored dynamic ranges from intrinsic brain pulsations at 1 Hz to ultralow-frequency vibrations at 5, 6.25, 7.8 and 10 Hz to the normal frequency range of MRE of 40 Hz. Surprisingly, we observed variations in brain stiffness over more than two orders of magnitude, suggesting that the in vivo human brain is superviscous on large scales with very low shear modulus of 42±13 Pa and relatively high viscosity of 6.6±0.3 Pa∙s according to the two-parameter solid model. Our data shed light on the crucial role of fluid compartments including blood vessels and cerebrospinal fluid (CSF) for whole brain properties and provide, for the first time, an explanation for the variability of the mechanical brain responses to manual palpation, local indentation, and high-dynamic tissue stimulation as used in elastography.

4D-CT Hyper-Elastography Using a Biomechanical Model.

Low dose computed tomography (LDCT) is the current gold-standard for lung cancer diagnosis. However, accuracy of diagnosis is limited by the radiologist's ability to discern cancerous from non-cancerous nodules. To assist with diagnoses, a 4D-CT lung elastography method is proposed to distinguish nodules based on tissue stiffness properties. The technique relies on a patient-specific inverse finite element (FE) model of the lung solved using an optimization algorithm. The FE model incorporates hyperelastic material properties for tumor and healthy regions and was deformed according to respiration physiology. The tumor hyperelastic parameters and trans-pulmonary pressure were estimated using an optimization algorithm that maximizes similarity between the actual and simulated tumor and lung image data. The proposed technique was evaluated using an in-silico study where the lung tumor elastic properties were assumed. Following that evaluation, the technique was applied to clinical 4D-CT data of two lung cancer patients. Results from the evaluation study show that the elastography technique recovered known tumor parameters with only 6% error. Tumor hyperelastic properties from the clinical data are also reported. Results from this proof of concept study demonstrate the ability to perform lung elastography with 4D-CT data alone. Advancements in the technique could lead to improved diagnoses and timely treatment of lung cancer.

A Tissue Mechanics Based Method to Improve Tissue Displacement Estimation in Ultrasound Elastography.

Cancer is known to induce significant structural changes to tissue. In most cancers, including breast cancer, such changes yield tissue stiffening. As such, imaging tissue stiffness can be used effectively for cancer diagnosis. One such imaging technique, ultrasound elastography, has emerged with the aim of providing a low-cost imaging modality for effective breast cancer diagnosis. In quasi-static breast ultrasound elastography, the breast is stimulated by ultrasound probe, leading to tissue deformation. The tissue displacement data can be estimated using a pair of acquired ultrasound radiofrequency (RF) data pertaining to pre- and post-deformation states. The data can then be used within a mathematical framework to construct an image of the tissue stiffness distribution. Ultrasound RF data is known to include significant noise which lead to corruption of estimated displacement fields, especially the lateral displacements. In this study, we propose a tissue mechanics-based method aiming at improving the quality of estimated displacement data. We applied the method to RF data acquired from a tissue-mimicking phantom. The results indicated that the method is effective in improving the quality of the displacement data.

Analytical Estimation of Out-of-plane Strain in Ultrasound Elastography to Improve Axial and Lateral Displacement Fields.

Many types of cancers are associated with changes in tissue mechanical properties. This has led to the development of elastography as a clinically viable method where tissue mechanical properties are mapped and visualized for cancer detection and staging. In quasi-static ultrasound elastography, a mechanical stimulation is applied to the tissue using ultrasound probe. Using ultrasound radiofrequency (RF) data acquired before and after the stimulation, the tissue displacement field can be estimated. Elasticity image reconstruction algorithms use this displacement data to generate images of the tissue elasticity properties. The accuracy of the generated elasticity images depends highly on the accuracy of the tissue displacement estimation. Tissue incompressibility can be used as a constraint to improve the estimation of axial and, more importantly, the lateral displacements in 2D ultrasound elastography. Especially in clinical applications, this requires accurate estimation of the out-of-plane strain. Here, we propose a method for providing an accurate estimate of the out-of-plane strain which is incorporated in the incompressibility equation to improve the axial and lateral displacements estimation before elastography image reconstruction. The method was validated using in silico and tissue mimicking phantom studies, leading to significant improvement in the estimated displacement.

A Composite Material Based Neural Network for Tissue Mechanical Properties Estimation Toward Stage Assessment of Infarction.

Cardiac biomechanical modelling is a promising new tool to be used in prognostic medicine and therapy planning for patients suffering from a variety of cardiovascular diseases and injuries. In order to have an accurate biomechanical model, personalized parameters to define loading, boundary conditions and mechanical properties are required. Achieving personalized modelling parameters often requires inverse optimization which is computationally expensive; hence techniques to reduce the multivariable complexity are in need. Presented in this paper is the fundamental blueprint to create a library of scar tissue mechanical properties to be used in modelling the healing mechanics of hearts that have suffered acute myocardial infarction. This library can be used to reduce the number of variables necessary to capture the scar tissue mechanical properties down to 1. This single parameter also carries information pertaining to staging of the scar tissue healing, predict its rate, and predict its collagen density. This information can be potentially used as valuable biomarkers to adjust existing or develop new treatment plans for patients.

Constitutive modeling of menisci tissue: a critical review of analytical and numerical approaches.

Menisci are fibrocartilaginous disks consisting of soft tissue with a complex biomechanical structure. They are critical determinants of the kinematics as well as the stability of the knee joint. Several studies have been carried out to formulate tissue mechanical behavior, leading to the development of a wide spectrum of constitutive laws. In addition to developing analytical tools, extensive numerical studies have been conducted on menisci modeling. This study reviews the developments of the most widely used continuum models of the meniscus mechanical properties in conjunction with emerging analytical and numerical models used to study the meniscus. The review presents relevant approaches and assumptions used to develop the models and includes discussions regarding strengths, weaknesses, and discrepancies involved in the presented models. The study presents a comprehensive coverage of relevant publications included in Compendex, EMBASE, MEDLINE, PubMed, ScienceDirect, Springer, and Scopus databases. This review aims at opening novel avenues for improving menisci modeling within the framework of constitutive modeling through highlighting the needs for further research directed toward determining key factors in gaining insight into the biomechanics of menisci which is crucial for the elaborate design of meniscal replacements.

Estimation of the hyperelastic parameters of fresh human oropharyngeal soft tissues using indentation testing.

Patient-specific finite element (FE) modeling of the upper airway is an effective tool for accurate assessment of obstructive sleep apnea (OSA) syndrome. It is also useful for planning minimally invasive surgical procedures under severe OSA conditions. A major requirement of FE modeling is having reliable data characterizing the biomechanical properties of the upper airway tissues, particularly oropharyngeal soft tissue. While some data characterizing this tissue's linear elastic regime is available, reliable data characterizing its hyperelasticity is scarce. The aim of the current study is to estimate the hyperelastic mechanical properties of the oropharyngeal soft tissues, including the palatine tonsil, soft palate, uvula, and tongue base. Fresh tissue specimens of human oropharyngeal tissue were acquired from 13 OSA patients who underwent standard surgical procedures. Indentation testing was performed on the specimens to obtain their force-displacement data. To determine the specimens' hyperelastic parameters using these data, an inverse FE framework was utilized. In this work, the hyperelastic parameters corresponding to the commonly used Yeoh and 2nd order Ogden models were obtained. Both models captured the experimental force-displacement data of the tissue specimens reasonably accurately with mean errors of 11.65% or smaller. This study has provided estimates of the hyperelastic parameters of all upper airway soft tissues using fresh human tissue specimens for the first time.

A Shape Optimization Technique to Predict Left Ventricle Ischemic Tissue Damage.

Damaged cardiac muscle tissue caused by ischemia leads to compromised cardiac function. While conventional imaging can view the ischemic tissue, currently there is no clinical way to quantitatively predict improved heart function after revascularization treatment. This increases the decision difficulty of treatment planning as there is no guarantee the heart function will improve enough to justify the cost of revascularization treatment. The complement of biomechanical modelling with conventional imaging offers an alternative method to determine the amount of ischemic tissue which can then be used as a potential predictor to estimate the range of functional improvement. A novel shape optimization technique is presented to predict the contractility of ischemic tissue in an in-silico left ventricle model that has suffered acute myocardial infarction. Preliminary results show that the proposed technique can reconstruct the damage caused by ischemic tissue within 18%. A range of minimum to maximum predicted cardiac improvement can then be given based on this error to help decide if the cost of revascularization treatment is justified.

4DCT Ventilation Map Construction Using Biomechanics-base Image Registration and Enhanced Air Segmentation.

Current lung radiation therapy (RT) treatment planning algorithms used in most centers assume homogeneous lung function. However, co-existing pulmonary dysfunctions present in many non-small cell lung cancer (NSCLC) patients, particularly smokers, cause regional variations in both perfusion and ventilation, leading to inhomogeneous lung function. An adaptive RT treatment planning that deliberately avoids highly functional lung regions can potentially reduce pulmonary toxicity and morbidity. The ventilation component of lung function can be measured using a variety of techniques. Recently, 4DCT ventilation imaging has emerged as a cost-effective and accessible method. Current 4DCT ventilation calculation methods, including the intensity-based and Jacobian models, suffer from inaccurate estimations of air volume distribution and unreliability of intensity-based image registration algorithms. In this study, we propose a novel method that utilizes a biomechanical model-based registration along with an accurate air segmentation algorithm to calculate 4DCT ventilation maps. The results show a successful development of ventilation maps using the proposed method.

Incorporating Pathology-Induced Heterogeneities in a Patient-Specific Biomechanical Model of the Lung for Accurate Tumor Motion Estimation.

Radiation therapy (RT) is an important component of treatment for lung cancer. However, the accuracy of this method can be affected by the complex respiratory motion/deformation of the target tumor during treatment. To improve the accuracy of RT, patient-specific biomechanical models of the lung have been proposed for estimating the tumor's respiratory motion/deformation. Chronic obstructive pulmonary disease (COPD) has a high incidence among lung cancer patients and is associated with heterogeneous destruction of lung parenchyma. This key heterogeneity element, however, has not been incorporated in lung biomechanical models developed in previous studies. In this work, we have developed a physiologically and patho-physiologically realistic lung biomechanical model that accounts for lung tissue heterogeneity. Four-dimensional computed tomography (4DCT) images were used to build a patient-specific finite element (FE) model of the lung. Image information was used to identify and incorporate inhomogeneities within the model. Mechanical properties of normal and diseased regions in the lung and the transpulmonary pressure driving the respiratory motion were estimated using an optimization algorithm that maximizes the similarity between the actual and simulated tumor and lung image data. Results from this proof of concept study on a lung cancer patient indicated improved accuracy of tumor motion estimation when COPD-induced lung tissue heterogeneities were incorporated in the model.

Estimation of the Young's moduli of fresh human oropharyngeal soft tissues using indentation testing.

Finite element (FE)-based biomechanical simulations of the upper airway are promising computational tools to study abnormal upper airway deformations under obstructive sleep apnea (OSA) conditions and to help guide minimally invasive surgical interventions in case of upper airway collapse. To this end, passive biomechanical properties of the upper airway tissues, especially oropharyngeal soft tissues, are indispensable. This research aimed at characterizing the linear elastic mechanical properties of the oropharyngeal soft tissues including palatine tonsil, soft palate, uvula, and tongue base. For this purpose, precise indentation experiments were conducted on freshly harvested human tissue samples accompanied by FE-based inversion schemes. To minimize the impact of the probable nonlinearities of the tested tissue samples, only the first quarter of the measured force-displacement data corresponding to the linear elastic regime was utilized in the FE-based inversion scheme to improve the accuracy of the tissue samples' Young's modulus calculations. Measured Young's moduli of the oropharyngeal soft tissues obtained in this study are presented. They include first estimates for palatine tonsil tissue samples while measured Young's moduli of other upper airway tissues were obtained for the first time using fresh human tissue samples.